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Abstract 


Maxwell's  equations  for  fields  with  sources  in  media  in  which  the  dielectric 
constant  and  permeability  are  unity  are  written  in  terms  of  a  spinor  notation 
which  resembles  the  one  used  for  Dirac's  equation  for  the  electron.  One  can  intro- 
duce Green's  functions  and  expansions  in  terms  of  complete  sets  of  orthogonal 
functions,  which  are  analagous  to  those  used  in  the  quantum  theory  of  the  electron, 
to  solve  Maxwell's  equations  in  more  compact  form  than  in  terms  of  the  conventional 
vector  notation. 

It  is  possible  to  express  both  real  and  complex  vector  potentials  simply  in 
tenns  of  the  field  strengths  by  means  of  this  notation.  Some  examples  will  be  given. 

Furthennore,  the  new  notation  enables  us  to  solve  in  a  simple  way  an  'inverse 
radiation  problem'  which  we  describe  as  follows: 

Consider  at  time  t  <  0  the  electromagnetic  field  to  be  zero.  At  time  t»0 
sources  are  turned  on  and  then  later  turned  off.  The  electrcwiagnetic  field »which 
results  after  this  process  has  been  completed^will  be  a  radiation  field.  We  can 
solve  the  problem  of  finding  the  nature  of  the  sources  which  will  lead  to  a 
prescribed  final  radiation  field.  It  is  shewn  that,  in  general,  the  sources  are  not 
unique  but  additional  conditions  can  be  given  which  will  make  them  so» 

In  Part  I  we  discuss  primarily  the  initial  value  problem  and  the  inverse 
problem  in  terms  of  field  strengths:  vector  potentials  play  no  role  whatever. 

In  Part  II  we  shall  introduce  various  vector  potentials.  As  is  well  known, 
it  is  necessaiy  to  use  such  potentials  when  one  wishes  to  discuss  the  interaction 
of  charged  particles  with  the  electromagnetic  field.  However,  in  contrast  to  the 
usual  procedure  in  which  one  obtains  differential  equations  for  the  vector 
potentials  and  then  obtains  field  strengths  in  terms  of  suitable  derivatives,  we 
shall  do  another  'inverse  problem':  we  shall  express  the  vector  potentials  in 
terms  of  the  field  strengths.  Vector  potentials  so  obtained  will  need  no  auxiliary 
conditions  such  as  the  Lorentz  condition  to  be  physically  meaningful. 
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1.     Introduction.  The  spinor  form  of  Maxwell's  equations 

Maxwell's  equations  in  free  space  with  sources  are: 


(1.1) 


curl  _E  +  ^  H  -  0, 
div  H  -  0  , 


div  E  "  Hnp  . 


(In   (1.1)  we  have  used  Gaussian  units  with  c  ■  1.) 

As  in  [ij   we  introduce  two  U-component  column  vectors  f  and  J? 


(1.2) 


(1.2a) 


T  - 


/M 


\^3 


^- 


iE, 


"2-^2 
H3  -  iE3 


(H^  -  H^,    H2  -  Hy,   H3 
and  )^l4  matrices  a       (i  «  0,1,2,3) 


lo 

I2 


H^,   etc.). 


P 
J. 


'3  ' 


0     T 

a  hI 


(1.3) 


'10  0  0 
0  10  0 
0  0  10 

.0  0  0  1, 


0     0-10 

0     0     0     i 

•10     0     0 

.0  -i     0     0 


0-100 
-10    0     0 
0    0     0  -i 
0     0     i     0, 


0     0     0-1 

0     0-10 

0     i     0     0 

•10     0     0 
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Maxwell's  equations  have  the  form 

(l.U)  -  1     y-  a^  ^  T  -  -  Unf 

where,   in  (l.U),   and  later 

(1.5)         X     »XQ«t,     X     --x^-x,     x^--X2"7,     T?'  - 


^  " 


z 


Equations  essentially  identical  with  (l.U)  have  also  been  given  in  [2] 
and  an  analagous  but  more  cumbersome  set  has  been  given  in  f3] .  In  \k\ 
and  pj  sets  of  equations  similar  to  (l.U)  have  been  given  for  the  truncated 
set  of  Maxwell's  equations  in  which  the  divergence  equations  have  been 
omitted. 

In  [l],  the  transformation  properties  of  the  vectors  T  and  J  under 
the  transformations  of  the  homogeneous  Lorentz  group  were  studied,  and  it 
was  shown  that  the  usual  transformation  properties  for  the  fields  and  sources 
are  obtained.  In  particular,  it  was  shown  that  the  field  function  T  transformed 
like  the  wave  function  for  a  spin  1  particle. 

In  the  present  paper  we  shall  study  the  solutions  of  (l.U)  and  their 
properties  in  terms  of  the  resemblance  of  (l.U)  to  Dirac's  relativistic 
equations  for  the  electron.  That  is,  we  shall  use  techniques  developed  for 
Dirac's  equations  to  obtain  solutions  of  Maxwell's  equations.  For  the  most 
part,  we  shall  obtain  familiar  results  in  more  compact  form.  However,  we 
shall  also  obtain  some  new  results  which  have  been  motivated  by  the  spinor 

notation.  For  example,  we  shall  be  able  to  find  all  time-dependent  sources 
which  lead  to  a  prescribed  final  radiation  field. 
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We  shall  now  develop  a  few  consequences  of  the  representation  (l.U).  By 
explicit  multiplication  we  note  that  the  operators  a  (i  ■  1,2,3)  satisfy  the 
following  multiplication  laws: 


(lo6) 


iah' 

-  ia'f 

-  ia')' 

aV. 

ia3-. 

ah"- 

aV- 

ia^.- 

a\^ 

a^al- 

ia^-- 

cV. 

Furthermore,  it  is  useful  to  note  that  these  operators  are  Hermitian. 

As  a  consequence  of  these  multiplication  rules  and  the  use  of  covariant 
and  contravariant  coordinates  as  given  by  (1.5),  one  obtains  the  following 
important  identity: 

^    Ha   \  /  T  3    1,   a  ^      /  T  j3    I,   ;i  \  /■  1  3 


This  property  is  very  similar  to  a  property  of  the  Dirac  operator.  On  applying 

1  ^   J  d 
-  xTL^"     S~"  *°  both  sides  of  (l.U)  we  obtain 

(1.8)      d'y.^Ec'^'^J. 

We  see  immediately  that  when  there  sure  no  source  terms,  the  vector  T  and 
hence  the  components  of  the  electromagnetic  field  satisfy  the  wave  equation. 
On  the  other  hand,  if  there  are  sources,  on  taking  the  top  component  of  the 
vectors  on  both  sides  of  the  equation.  (1.3)  and  using  T^^  5  0  ,  we  obtain  as 
a  necessary  condition  for  the  solution  of  (l.U),  the  equation  of  continuity: 


-  h  - 


''•''  tS'^^*'^^-'- 


In  most  of  our  applications  it  will  be  useful  to  separate  the  time 
from  the  space  variables  and  in  place  of  (l.U)  we  shall  write 

where  H-^  is  given  by 

The  operator  IL.  is  analagous  to  Dirac's  operator  for  the  kinetic  energy 
of  the  electron.  We  can  now  state  our  primary  objective:  we  propose  to  solve 
(1.10)  by  word-ng  essentially  in  the  spectral  representation  of  the  operator 
Hq  instead  of  using  the  usual  Fourier  transformations  in  terms  of  wave  numbers. 
Some  use  of  this  approach  has  been  made  in  (TJ  and  [5]  but  the  treatment  in 
these  references  is  incomplete  and  also  somewhat  cvunbersorae  because  the  equations 
analagous  to  (1»10),  which  these  authors  use, do  not  represent  the  complete  set 
of  Maxwell's  equations. 

Finally,  in  this  section  we  shall  define  suitable  inner  products  of 
column  vectors  and  show  how  the  energy  conservation  laws  may  be  derived  simply 
in  terras  of  these  definitions. 

Let  us  define  two  different  types  of  inner  products  of  column  vectors. 
Consider  two  vectors  A(x)  and  B(x)  which  are  given  by 


(1.12)   A(x)  =  \yr\  B(x) 


B^(x) 
B2(x) 

B3(x) 
In  (1.12),  and  later,  x  represents  collectively  the  three  dimensional  space 
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12   3 
coordinates  (x  ,  x  ,  x)»     The  first  type  of  inner  product  which  we  wish  to  define 

is  essentially  the  Hermitian  inner  product  in  finite  dimensional  vector  space. 

We  denote  it  by  A(x)«B(x) 

3    * 
(1.13)      A(x)-B(x)  '  y    A,  (x)  B, (x)  , 

1"0 
Where  the  asterisk  means  the  complex  conjugate.  The  second  inner  product,  which 
we  denote  by  (A,B),  is  the  Hermitian  inner  product  used  in  Hilbert  space 

3 


3   f  r 

(I.IU)   (A,B)  -  r;    A*(x)  B, (x)  dx  -   A(x).B(x)  dx  . 

i-o  J  J 


It  is  important  to  note  that  the  operators  a  are  Hermitian  with  respect 
to  both  inner  products,  i.e.. 


(1.15) 


A(x)»o^B(x)  -  a^A(x)»B(x)  , 
(A,a^B)  -  (a^A,B)  . 


We  shall  now  show  how  the  energy  density  and  Poynting  vector  may  be 
expressed  in  terms  of  these  inner  products  by  deriving  the  appropriate 
conservation  law.  Let  us  define  the  vectors  A(x)  and  B(x)  by 


(1.16) 


A(x) 


m 


-1 


T(x)  , 


j  a 


r-' "  ^  k "  i? ' 


and  take  the  inner  product  (A,B).  Because  of  Maxwell's  equations  (1»U)  we 
have 

(1.17)   (A,B)  -  ^  H  I  I(x)  «  a^   -^  T(x)  dx  -  |  J  Y(x)  •  J  (x)dx  . 


On  the  other  hand,  on  taking  the  product  (B,A),  we  see  that 
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(1.18) 


^fcf"'^'<~'-"~'^^ 


1 


$  (x)  •  T(x)  dx 


Since  the  operators  a  are  Hermitian ,  equation  (1.18)  becomes 


-^  T(x)  •  a^f  (x)  dx  -  -  i   5  (x)  .  T(x)  dx  . 


On  adding  (1.17)  and  (l.l8a)  we  obtain 


(1.19)   r 

where 


Xs^(x)  dx  ^~\   S°(x)  dx  * 
dx 


et 


3 


-4r  s^(x)  dx 
dx^   ^  " 


P(x)  dx  , 


(i  -  0,1,2,3) 


(1.19a)    I 


S^(x)  -  T(x)  •  a^f (x) 

3 
P(x)  -  Re  if(x)  •  f  (x)  -  7"  E,j. 

where  Re  means  "  real  part." 

It  is  easy  to  see  that  P(x)  is  the  time  rate  of  change  at  which  the  electric 

field  does  work  per  unit  volume  on  the  sources.  Hence  S  is  the  energy  density 

and  the  components  S  (i  ■  1,2,3)  are  the  components  of  the  Foynting  vector. 

Since  the  field  f  is  arbitrary  we  can  strip  of  the  integrals  in  (1.19)  and 

obtain  the  familiar  differential  form  of  the  conservation  of  energy 

,0 


(1.19b) 


3S 


i^  +  div  S  -  E»J 


Various  other  conservation  laws  can  be  obtained  in  an  analagous  way. 

2.  The 'eigenf unctions*  of  H-.;  the  x  representation;  the  P  representation 

It  is  our  objective  to  solve  the  initial  value  problem  of  Maxwell's 
equations  (l.li)  in  the  infinite  domain  using  techniques  suggested  by  the  resemblance 
of  (l.h)  to  Dirac's  equation  for  the  electron.  Accordingly,  we  shall  work 
in  the  spectral  representation  of  the  operator  IL..  That  is,  we  shall  expand 
the  column  vectors  I  and  J  in  terms  of  the  eigenfunctions  of  H^.  In  terras  of 
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this  expansion  Maxwell's  equations  (l.Ii)  are  solved  simply. 

Actually,  there  is  a  slight  complicaticn  in  the  program  because  of  the 
fact  that  the  top  component  of  Y  vanishes.  Hence,  for  the  purposes  of  expanding 
T,  we  should  like  to  distinguish  the  eigenfunctions  of  H^,  for  which  the  top 
components  vanish,  frcati  the  others.  However,  it  turns  out  that  there  are 
only  two  such  eigenfunctions  and  these  do  not  form  a  complete  set  which  span 
the  space  for  which  the  top  component  of  the  column  vector  is  zero.  Accordiniily, 
we  introduce  another  'eigenf unction'  which  is  orthogonal  to  the  acceptable 
e'.genfimctions  of  H-.  such  that  the  top  component  is  zero.  Then  all  field 
functions  Y  can  be  expanded  in  terms  of  these  three  'eigenfunctions. '  Finally, 
since  we  want  to  be  able  to  expand  vectors  J  which  have  top  components  also, 
we  introduce  a  fourth  'eigenf unction'  orthogonal  to  the  other  three.  The 
last  two  'eigenfunctions'  are  not  really  eigenfunctions  of  H_  -  only  the  first 
two  are  -  but  as  indicated  above  they  caaie   in  naturally  and  are  actually  related 
closely  to  eigenfunctions  of  H-  independent  of  the  acceptable  ones. 

Rather  than  give  a  detailed  derivation  of  the  'eigenfunctions'  of 
Hq,  we  shall  instead  give  their  forms  and  properties. 

We  shall  designate  the  four  'eigenfunctions'  by  X  (x|p, e)  where  p  is 
a  three-dimensional  vector  and  where  e  has  the  values  ♦  1,0,1'.  That  is,  the 
four  values  of  e  label  the  fotxr  'eigenfunctions. '  The  eigenfunctions  are 
four  component  column  vectors  which  we  may  also  write 


(2.1)  X(x|p,e) 


;Lo(x|p,e) 


Explicitly,  the  set  of  'eigenfunctions'  which  we  use  is 
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1£.X 


7(x|p,e)  -  y75 


-  (•?  ^  -  ie*?  ) 
'  X  'z      y 


<^"'"'  vf-e^ 


-  {-)  y  ^,  ♦  ie1  ^) 


for  e  ■  ±1, 


2    2 


(2.2) 


< 


ip'X, 
9(.(x|p,0)  - Cl«. 


(2n)' 


V. 


In  (2.2),  "^f  is  the  unit  vector  in  the  direction  of  p,i.e., 
(2.2a) 


P 


where  p  »  Ipj.  Also  £«j  is  the  usual  three-dimensional  scalar  product  of  p  and  x. 

It  is  easilv  verified  that  the  'eigenfunctions •  are  orthonormal  to  each 
other,  i.e.,  they  satisfy 

(2.3)  X-     P^i^^'P*^^  \(x|£',c')dx  -  j  X(x|p,c).X(x|p',e')dx  -  6(p-p')6^^,  , 


where 


(2.3a) 


6   ,  -  0  , 
e,  e'     * 

6    •  1  . 


This  set  of  'eigenf unctions'  also  satisfies  the  completeness  relation: 
(2.1.)    YL   Ui(x!p,e)  X,(x'|p,e)dp  -  6(x-x06,, 

where  ^  means  that  the  summation  is  taken  over  all  the  ' eigenf uncti ons. ' 
e 
As  a  consequence  we  can  expand  any  four-component  column  vector  A(x) 
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A(x) 


in  the  following  way 


(2.?)  A(x)  -  Yi 
e  J 


X(x|p,e)  A(p,c)dp 


or 


where 


AiW  -e)  Xi^b^^  A(p,e)dg, 


(2.6)  A(^,e)  -F    X^(x|p,e)  A^(x)dj  -   X  (x|p,  f.)»A(x)dbc  e 

If  we  choose,  we  may  regard  A. (x)  and  A(p, c)  as  representing  an  abstract 
vector  A  in  Hilbert  space  in  two  different  representations,  which  we  may  call 
the  x-representation  and  the  g-representation, respectively.  The  relations 
(2.^)  and  (2.6)  constitute  a  transfonnation  between  the  two  representations.  It 
is  easy  to  show  that  the  transformation  preserves  the  inner  product  in  the 
Hilbert  space  sense.  That  is,  let  A. (x)  and  A(p, e)  represent  the  sane  vector 
A  in  the  x-  and  p-representations,  respectively.  Likewise,  let  B. (x)  and 
B(p, e)  represent  the  abstract  vector  B  in  the  two  representations.  Then,  from 
(2.5)  and  (2.6)  together  with  the  completeness  and  orthononnality  relations 
(2.3)  and  (2.U),  one  can  show 


3   f  » 

(2.7)   (A,B)  -  5~    A. (x)  B. (x)dx  -  J~ 

iTo  J  6  J 


A  (pjc)  B(p, 6)dp 


The  transformation  functions i (x I p,il)  are  true  eigenfunctions  of  H„ 
which  have  their  top  components  zero.  They  satisfy 

(2.8)   fL,X(x|p,e)  -  epX(x|p,£)  ,   (e-=±l)  ,     where  p  -  |p|   . 


-  10  - 


The  'eigenf unction'  llCxIpjC)  is  not  a  tn^e  eigenf unction  of  H^.  However, 
as  shown  below,  H^  acts  very  simply  on  it.  This  is  the  transformation  vector 
which  is  introduced  to  enable  us  to  expand  all  electromagnetic  field  functions 
7. (x)  such  that  I^(x)  -  0. 

(2.9)  H-.  l(x|p,0)  -  -  p  l(x|p,  ?)  ,     (p-IpI). 

The  transformation  function  7'(x|p,  ?■)  is  introduced  to  enable  us  to  ex- 
pand any  four-component  colvnnn  vector.  The  operator  Hq  acts  in  a  simple  fashion 
on  this  vector  also. 

(2.10)  Hn  l(x|p,  ?)  -  -  p   X(x|p,0). 

Another  relation  that  will  prove  very  useful  is  the  divergence  property 
r 


(2.11) 


r  A  ^i(sip»i^)-o  , 

1=1  dx 

f  A  A,(x|p,o).^£57j    e^E-^ 

I?L  d?  ^  ^  ^              (2n)^/2 


As  we  shall  show,  the  relations   (2.11)  enable  us  to  separate  the 
longitudinal  from  the  transverse  field  in  a  simple  fashion. 

We  are  now  able  to  expand  otir  time-dependent  electromagnetic  field 
vector  f(x:t)  and  sources  $(xjt)  in  terms  of  the   'eigenvectors'    7-   (x|p, e) 


(2.12) 


where 


(2.12a) 


T(x;t)  -  YL 


e    J 


< 


$(xjt)  'YL 
IT  J 


%(xjp,e)  f(p,ejt)dp  , 

t(xlp,f)  J(p,ejt)dp  , 
"  ^  ^,  '-' 


I(p,e|t)  -P       XiU|p,e)  T^(x|t)dx  , 
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Since  we  require  T-^(x;t)  =  0,   we  see  from   (2.2)  that  in  the  £7 
representation  we  must  have 

(2.13)  lipftit]  a  0   . 

Hence  we  may  write 

(2.1h)         T(x{t)  -  T'^(x;t)  +  Y^(x;t)   , 

/ 

X(x|p,e)  T(p,e|t)dp  , 


v^ere 


e=±l 


'V 


■^' 


I^(xit)  -       X(x|p,0)  T(p,0;t)dp  . 


Now  Y     is  the  transverse  field,   since  from   (2.11)  we  have  J~    — ^  T.    "  0  , 

T  T  L  ^"^  ^* 

which  leads  to  div  E     ■  div  H     ■  0.     Likewise  T     is  the  longitudinal 

part  of  the  field.     From  the  orthogonality  relations   (2.3)  we  see  that  the 

longitudinal  and  transverse  fields  are  othogonal   (in  the  Hilbert  space  sense)  to 

each  other: 

(2.15) 

Because  of  (2.15)  the  energy  of  the  field  splits  up  into  a  longitudinal 
energy  and  a  transverse  energy:  The  energy  En  over  all  space  is  given  by 


I  Y^(xjt).Y^(3yt)dx  -  0  . 


(2.16)     En 


+  ^      T^(xjt)«T^(xjt)dx  "  Erl  +  E?l  . 


We  now  consider  the  sources.     The  equation  of  continuity  which  is  a 
necessary  condition  for  the  solution  of  Maxwell's  equations 


(2.17) 


3      a$ 
i^o  dx^ 


lead  to,   on  using   (2.11)  and  (2.12),   a  restriction  on  J(p, e;t),   namely 


-  1?  - 


dj(p,r,t) 

(2.1f;)         J(p,0;t)  -  i    —^ . 

Hence  in  the  p-representation  (2.l8)  replaces  the  equation  of  continuity. 
The  condition  that  J. (xit)  is  real  leads  to  necessary  and  sufficient 
s-vmmetry  conditions  on  ^(p,  e;t),   namely 

^■^(P^tljt)     »    |(-     P,±ljt)     , 

i*(p,T^it)  =  j(-  g,rit) , 

f*(p,0;t)     -  -  $(-  PiOjt)   . 


(2.19)  < 


We  shall  find  (2.19)  of  especial  importance  in  the  'inverse  problem*  of 
determining  sources  from  radiation  patterns. 

3.  The  solution  of  the  initial  value  problem  without  sources;  circularly  and 
plane  polarized  waves j  the  Green's  function  for  the  initial  value  problem; 
energy  density  in  the  jaj-representation 

We  shall  now  solve  the  initial  value  problem  for  Maxwell's  equations 
without  sources,  using  expansions  in  terras  of  the  'eigenfunctions '  A-(x|p, e). 
Later  we  shaii  extend  our  results  to  the  cases  where  sources  exist. 
Maxwell's  equations  are: 


(3 


'1^     (-T  it-''^o)'-° 


As  in  S3Ction  2  we  write 


(3.2)        Y(x;t)  -7"    j    %(x|p,c)  f(p,e;t)dp 

anri  siibstitute  into   (3.1).     Equations   (2.8)   -  (2.10)  shovr  u£  how  Hq  acts 
on    X (x I p, e )  and  we  obtain 
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(3.3) 


-tH   i(x|p,e)  |r  (p,e;t)dp-  XI    ^P  /t(x|p,e)  Y(p,e;t)dp 


pX(x|p,  ?-)  T(p,Ojt)dp  + 


pX(x|p,0)  Y(p,?;t)dp  =  0  . 


Since  the  'eigenf unctions •  are  orthonormal  we  may  set  the  coefficients  of 
X(x|p,  e)  equal  to  zero  and  obtain 


^  a¥(p,e;t) 
(3.'i)    -  J  '^ epY(p,6;t)  -  0  , 


(e  »i  1)  , 


(3.'4a) 


aT(p,0;t) 
i— -^ *  pl(p,?"it)  -  0 


I  — ft *  PT(p,0;t)  -  0  , 


As  mentioned  in  Section  2,  the  condition  that  T-.(xjtj  S  0  leads  to 

(3.5)  Y(p,?';t)  5  0 

and  hence  from  the  second  equation  of  (3,Ua), 

(3.6)  T(p,0;t)  £  0  . 

As  a  consequence  of  (3.6 j,  we  see  that  the  solutions  of  Maxwell's 
equations  without  sources  lead  to  transverse  fields  only. 

The  equations  (3.U)  which  give  the  transverse  fields  are  easily  solved. 
The  solutions  are 

(3.7)  I(p,e,t)  -  e'^'^(^'^0)   Y(p,e;to)  ^   (^  „.^  ^)  ^ 

where  Y(p, e|t-)  is  the  value  of  Y(p, e;t)  at  time  t=t^.  Hence  the  solution 


of  Maxwell's  equations  without  sources  is 

(         -iep(t-tQ) 
(3.8)  Y(x;t)-E   U(x|p,e)e         ^^(p,  e,-t^)dp 


lu  - 


when  we  give  the  value  of  T(x;t)  at  t  «  t-^  in  the  ^-representation.  We  shall 
shortly  show  how  ¥(x;t)  may  be  obtained  directly  in  terms  of  Y(xit_)  through 
the  use  of  an  appropriate  Green's  function. 

However,  we  shall  first  show  that  the  solution  (3.9)  consists  of  a 
superposition  of  circularly  polarized  waves.  That  is  to  say,  the  solution 

(3.8)  consists  of  a  superposition  of  solutions  Y(xjt|p,£)  where 

(3.9)  Y(xjt!p,e)  -  t(x!p,e)e'-^^P*  ,  e  -  4  1  . 

We  maintain  that  the  solutions  T(xjt|p, e)  represent  circularly  polarized 
electromagnetic  waves  with  frequency  p  =  |p(  travelling  in  the  direction  p 
for  e  "  1  and  travelling  in  the  direction  -  p  for  e  -  -  1.  To  prove  this 
assertion  let  us  choose  our  coordinate  system  such  that  the  x-axis  coincides 
with  the  vector  p.  Then^on  using  the  explicit  forms  for^(x|p,e)  given  by 
(2.2)  and  noting  that  ^      -  1,  7   -  •/  -  0,  we  have 


(3ao,  „x,t,p,o  -  ^  ^  .^^<--"  ° 


-JLfc 


Since  the  electric  field  strength  is  the  negative  of  the  imaginary  part  of 
Y,  while  the  magnetic  field  strength  is  the  real  part,  we  have 

E  "  0  H  -  0 

X  x 

(3.11)   E  -  ^  A  lA  cos  p(x  -  et)       H  -  "  ^  ^yA  sin  p(x  -  ct) 

E  » i«-7».  sin  p(x  -  et)      H ^75-  cos  p(x  -  et)  • 

V/?  (2n)3/2  z   ^  (2n)^/^ 

These  are  just  the  usual  forms  for  circularly  polarized  eiectrcwaagnetic  waves 
which  travel  in  the  positive  x  direction  for  e  ■  1  and  the  negative  x  direction 
for  e  ■  -  1. 
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It  is  interesting  to  note  that  our  representation  leads  naturaxly  to 
a  decomposition  in  terms  of  a  sum  of  circularly  polarized  waves  rather  than 
plane  polarized  waves.  In  fh]   and  f^   a  similar  decomposition  is  used. 

One  can,  of  course, construct  plane  polarized  waves  from  circularly 

polarized  waves.  Consider  the  solution  of  Maxwell's  equations  f '(xjt|p, e) 

defined  by 

Y(x;t!p,f.)  +  e¥(x,-t|-  p,-  e) 
(3.12)   Y'(x;t|p,6)  -  — ^-^^. ^-^-^ :^- 

We  m? intain  that  T'  represents  a  plane  polarized  wave  with  frequency  p  ■  |p| 
traveling  in  the  direction  p  for  e  ■  1  and  in  the  direction  -  p  for  g  =  -  lo 
Let  us  again  choose  the  coordinate  system  so  that  the  x  axis  coincides  with 
the  vector  p.  Then  by  taking  the  real  and  imaginary  parts  of  f '  we  obtain 

E     =  0  H     «  0 

X  X 

\r ^     o/o  -^^  cos  p(x-et)     H^  -  -»^k-^5^  i^L^l  sin  p(x-et) 

y      sn  i2ny'  ^      .y?  (2n)^/2       ^ 


z 


1  (e-l)      •        /       4.\     T,  1  (e+1)  /       .X 

— T-TK-  -Hr^  sm  p(x-ct)     H     ■=  -^ — 77^  >.-..,    cos  p(x-et) 

•y?  (2n)3/2       2  z     j^  (2n)3/^  ^ 


from  which  our  assertion  is  proved. 

We  shall  now  express  the  general  solution  of  T(xjt)  of  Maxwell's 
equations  in  terms  of  the  initial  value  I(xit  )  through  the  use  of  an 
appropriate  Green's  function.  We  may  write  the  solution  (3.8)  as 

(3.11;)  T,(x|t)-E   I  X.(x|p,e)e~^^P^*"*0^  T(p,e|t^)dp. 
^  "      e»ll  J   1  A(  ^  ,^         ^     ~^ 

According  to  (2.32a),  W3  have 
(3.15)    Y(p,e,-tQ)  -  r  I  X*(x'|p,e)  Yj(x';tQ)dx'. 
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On  substituting  (3.l5)  into  (3. Hi)  we  obtain 


(3.16)  I.(xjt)  «  L, 


Xi(x|p,e)/.;(x.!p,Oe"^'''  """'^  ^j^^S'^S^^'^i^ 


which  we  may  also  write  in  terms  of  a  Green's  function  notation  as 
(3.16a)  T(xjt)  -   G(x|t|x';tQ)  ¥(x';tQ)cbc' 

where  G(xjtlx'jt')  is  an  operator  Green's  function  whose  elements  are  given  by 

(3.17)  G,,(x;t|x';t')  -  r     X.(x|p,Ol!(x'ip,e)  e"''^^*"*'W  • 

By  using  the  explicit  forms  of  the  eigenf unctions  which  occur  in 
(3.17)  as  given  by  (2.2)  we  can  evaluate  the  matrix  elements  of  the  Green's 


function  G.  We  find 


fooj^i**!^''*')  -0     ' 


(3.13) 


G  (x,t|x',t')  - 
^G^2(^t|x',t')  - 
and  cyclically 


2        ^  i  2 


dt    dx 
-  1 


ata?  dxdx 


(j-0,1,2,3) 

B(x|tix';t')  ,   (i-1,2,3) 

B(x;t!x','t') 


*^     /^* 


G  (x,t|x',t')  -  G*  (x,t|x',t') 


where 


(3.lBa)    B(x;tlx';t') 


fdx-x'i^.  it-t-r) 


The  function  (x)  is  the  Heaviside  step  function  defined  by 


(3.18b)    '/(x)-! 
-  0 


X  >  0 
X  <  0 
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The  first  of  equations  (3.18)  follows  from  the  fact  that  the  top 
cwnponent  of  the  eigenf unctions  involved  in  calculating  G  is  aero,  We  shall 
prove  the  second  equation  of  (3.18)  for  i  ■  1.  The  remaining  equations  are 
proved  similarly. 

Frcm  (3.17)  and  (2.2)  we  have 

(3.19)  G,,  (xjt[x';f)  -  Z         /i(^!P.e)  X*(x|p,e)e-^^P^*-^'W 


2(2n)3  J 


"(20^  J 


X   'y 

•(x-x') 

/  X  /y 


'  X   'y 

r    ip«(x-x') 

e  ^ 


-JL-.   e-'^'^^'^'^U  -O^lcos  p  (t-t')dp 
(2„)3  J         L   '  xJ 


-  3'  +  ^^ 

3?   3x  /  (2n)' 


Let  us  now  evaluate  the  integral.  Let  us  set 


ip«(x-x') 

»  rv   A^  -v 


cos  p(t-t')dp  • 


(3.20)   ?  -  X  -  X'  ,  C  -  kl  ,  X  -  t  -  t'  . 

V     ~     '^  'V' 

Then,  on  using  polar  coordinates  with  x-x'  coinciding  with  the  p  axis 


we  have 


(3.21) 


[  ip-(x-x') 
1    e  ^  ^^   ^^g  p(t-t')dp 


2n         00 


(2n) 


(2n) 


i-^  J  d0  I  dp  cos  pX  j  sin  ©  e^P^  °°^  ^de 


0    0 


00 

2   f   d£ 


A     — '^  cos  pX  sin  p4 
(2n74  i   P 


(2n)^4 


+  00 


cos  pX  sin  p5  • 


-00 
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This  integral  can  easily  be  evaluated  by  the  usual  methods. 
One  shows  that  if  X  >  0 


1 

p^  cos  pX  sin  PC  >  jj^  , 

00 

4  >  X 

(2n)2^  J 

-  0   , 

C  <  >^  J 

while,  if  X  <  0, 

,+  00 


T. 


d£ 


cos  pX  sin  p4  " 


(2n)'4  j.,   ^    —  —   115^  ' 


•00 


0 


4  >  -X 
4  <  -X  . 


These  results  can  be  summarized  by 


(3.22) 


(2n)2? 


.TUU 

dp 
P 

cos 

pX 

sin 

P4 

.  1 

7(4- 

M) 

-00 

-  1 

f(«'- 

x2) 

-  3(x,t|x',t')  . 

To  conclude  this  section  we  shall  discuss  the  ^-representation. 

For  Maxwell's  equations  without  sources  the  p-representation  has  a 
particularly  simple  interpretation,  since,  as  a  consequence  of  the  length-preserving 
character  of  the  transformation  between  the  x-  and  p-representations,  the 
energy  of  the  transverse  field  may  be  written  as 


(3.23)   En  «  ^ 

e=±l  , 


T  (p,ejt)  Y(p,e;t)dp 


In  the  present  case  we  have  seen  that  the  solutionsof  Majcwell's  equations 
are  a  superposition  of  circularly  polarized  waves  in  various  directions.  Hence, 
from  (3.23)  we  see  that  T  (p, e;t)  T(p, ejt)  can  be  interpreted  as  the  energy 
density  in  p  space  of  circularly  polarized  radiation  oT  frequency  p  "■  |p!  prop- 
agated  in  the  direction  ep. 
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U.   The  solution  of  the  initial  value  problem  with  sources.  Separation  of 
the  longitudinal  and  transverse  fields* 

We  shall  now  use  the  expansion  of  f  and  ^  in  terms  of  the  complete 
set  of  'eigenf unctions'  to  solve  Maxwell's  equations  with  sources.  As  in  Section 
2,  let  us  write 


(U.l) 


T(x;t)  =  J2 
e 

$(x,.t)  =  j:: 


e  J 


"XCxIpje)  T(p,f.;t)dp  , 
/((x|p,e)  f(p,  e;t)dp  . 


On  substituting  Y(xjt)  and  J(x,t)  as  given  by  (U.l)  into  Maxwell's 
equations 


(^•2)     [-^|p-Ho]T(x;t)-J(x;t)  , 

and  using    (2.8)-(2.10),   we  obtain 

aY(p,e;t) 


(U.3)   -In 


l(x|p,£)  g^   '         +    I    X(x|p,0)   pT(p,  r;t)d£ 

■ 

■X(x|p,t)pT(p,0;t)dp-  7~       ep?f(x|p,  e)Y(p,  eit)dp 
~  ^"  H      1(^1  P,e)  |(p,e;t)dj)  . 


By  identifying  coefficients  of  the  ' eigenf unctions '  X(x|p, e)  we  obtain  equations 
for  T(p,e;t)  in  terms  of  f(p,ejt),  i.e.,  for  T  in  terms  of  f  in  the  p-representation. 
The  equations  for  the  transverse  part  of  the  field  are  obtained  by  identifying 
the  coefficients  of  X(x|p,±l)  and  we  find 

(b.U)  "  J  5^  T(p,f.it)  -  ep(p,fjt)  "  -  Un|(p,ejt)  ,   (e  -  +1)  , 

As  before,  we  require 
(U.5)  Y(p,  ^;t)  m   0  . 
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Hence,  by  identifying  the  coefficients  of  X  ixl^,?^) 

(h.6)    Y(p,o;t)  -  -  kn   |(p,r,-t)  , 

and  by  identifying  the  coefficients  of  X  (xIp,C) 

(U.7)  -  ^g|  (p,0;t)  -  -  Un|(p,0;t)  . 

In  order  that  I(p,C;t)  be  given  both  by  (U.6)  and  (U,7)  we  obtain  a 
relation  that  J  must  satisfy: 

sjCpj-r  jt) 

which  is  just  equation  (2.18):  this  was  seen  to  be  equivalent  to  the  equation  of 
continuity 

i»0  ex 

In  our  notation  (as  in  the  usual  treatment)  the  transverse  and  longitudinal 
fields  uncouple.  The  transverse  field  in  the  p-representation  is  given  by  (U.U). 
Because  of  the  first  derivatives  in  time,  one  can  consider  an  initial  value 
problem  which  this  poses. 

On  the  other  hand,  the  longitudinal  field  itself,  rather  than  its 
derivative,  depends  on  the  sources  and  is  therefore  simpler  to  obtain  than  the 
transverse  field.  No  initial  value  problem  is  involved. 

Since  the  relation  (h.fc)  for  the  longitudinal  field  is  so  simple,  we 
shall  solve  it  in  terms  of  the  x-representation. 

On  using  the  expressions  for  the  eigenf unction  (2.2), 
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(li.lO)  I^Uit)  -  -  Un    ^  p"^ l(p,  ?-;t)d^ 


(2 


nip<x-x • )     . 
S-^^i^  <iP' f  iJo<J''"*S'   .    i-  1.2.3, 


where  ^  •''1    >  etc.  Sine©  $Q(x;t)  »  p(x;t)  we  see  that 

(U.ll)   n(x|t) 5-  T      ; 

J  '"       2n^  ax-'  J  J     p' 
But 


dp  p(x';t)dx« 


1 


— ? ^B  ■  — 


x-x,' 


Hence 


p(xjt) 


x-x« 


dx»  ,    i  -  1,2,3  . 


(U.12)  Y.(xjt)  -  i  -^ 

Since  f  J"  »  hJ^  -  iE^  , 

we  obtain  the  familiar  result  that  the  longitudinal  magnetic  field  is  zero. 

(U.13)  H^  2  0, 

whereas  the  longitudinal  electric  field  is  related  to  charge  density  in  the  same 
way  as  the  electrostatic  field: 


(U.lU)   E(x;t)  -  -  grad 


1 


p(x';t) 


dx'  . 


For  the  sake  of  completeness  we  shall  also  write  the  energy  of  the 
longitudinal  waves 

(U.15)  Eft  -  ^  r  J  Y5;*(x|t)  ^\U}t)d^  -  l^  J  T*(p,Ojt)  T(p,Oit)d^  , 
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which  follows  from  the  length-preserving  character  of  the  transformation  from  the 
X  to  p  transformation.  On  using  (U.6)  and  (2.6)  and  finally  (U.12),  then 


(U.36)  En  -  2n  I  f*{p,T}t)   f(p,7-;t)dp 

^   dx  p(xit)   dx'p(x';t) 

^  r  r  dx  dx'  p(x;t)  p(x«;t) 
"^J  J      Ix-x'l 


_2i 
(21 


i^» (x-x' ) 
dp£ ^ 

~    P 


which  is  the  usual  expression  for  the  energy  of  the  longitudinal  field. 

Having  disposed  of  the  longitudinal  field,  we  shall  now  discuss  the 

initial  value  problem  for  the  transverse  field. 

The  general  solution  of  the  differential  equation  (U.U)  is 

(U.17)   T(p,e;t)  -  e'^^P^^'^ohip,t}tQ)  *   Unie'^'P^    J(p,ejt  •  )e^'P*'dt,   (e«±l)  . 

^0  " 
The  first  term  on  the  right  represents  a  solution  of  Maxwell's  equation  without 
sources,  while  the  second  term  shows  the  effect  of  sources.  It  should  be 
noted  that  $^(x',t)  -  p(x;t)  has  no  effect  whatsoever  on  the  transverse  field. 

We  can  also  write  the  solution  (U.l?)  in  terms  of  the  x-representation. 
We  have, as  a  consequence  of  (2.12)  and  (2.12a), 

Tj(x,'t)  -21     Xi(?,IP»e)  f'^(p,e,-t)dp 
(U.18)        -  Y:       f  J  Xi(x|p,e)X*(x'|p|e)e-^^P(*-*O^Tj(x',tQ)dx'd£ 

*  Uni  j  dt'  II^I   Xi(x|p,e)Xj(x'|p,e)e-^^P(*-*'^|j(x',t')dx'd£. 


Hence 
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(U.19)  Y'^Ujt)  -  ]  djc'G(xjt|^',-tQ)  T^(x';t^)  +  Uni    dt'  I  cte'Q(xjt [x'jt • )  J(x«,-t') 

where  G(xjt|x'jt')  is  the  Green's  function  given  by  (3.17)  and  (3.18).  Equations 

T  L 

(U.19)  for  T  and  (U.12)  for  Y  represent  the  complete  solution  of  the  initial 

value  problem  for  Maxwell's  equations  with  sotirces. 

5»   The  influence  on  radiation  of  sources  which  are  switched  on  for  a  finite 

time.  The  problem  of  finding  sources  which  will  give  a  prescribed  radiation 
field. 

In  the  previous  section  we  have  solved  the  problem  of  obtaining  the 
electromagnetic  field  in  terms  of  the  sources  in  the  time-dependent  formulation. 
In  the  present  section  we  want  to  consider  a  special  situation  problem,  namely 
the  case  in  which  the  sources  are  identically  zero  except  for  a  finite  time  interval. 
Such  a  situation  leads  to  a  special  initial  value  problem. 

Before  the  sources  are  switched  on,  the  only  solutions  of  Maxwell's 
equations  are  those  given  in  Section  3  irtiich  correspond  to  electromagnetic  radiation. 
In  the  interval  during  which  the  sources  are  on,  one  obtains  a  more  complicated 
transverse  field  and  also  a  longtitudinal  field.  After  the  sources  are  switched 
off,  one  again  obtains  solutions  corresponding  to  electromagnetic  radiation,  which, 
however,  differ  in  general  from  the  initial  fields  before  the  sources  were 
switched  on. 

A  physically  interesting  problem  is  to  obtain  the  final  field  from  the 
initial  field  and  sources.  The  problem  is  analagous  to  the  scattering  problem 
in  quantum  mechanics.  Since  we  have  given  the  general  solution  of  Maxwell's 
equations  in  the  previous  section,  we  can  easily  obtain  the  final  field  from  the 
initial  field  and  the  sources. 

We  shall  also  be  able  to  solve  the  'inverse'  problem  which  may  be 
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described  in  the  following  way:  we  prescribe  the  initial  field  and  the  final 

field.  We  are  required  to  find  the  sources  which  lead  from  the  initial  field 

to  the  final  field.  A  particular  case  is  that  which  occurs  when  the  initial 

field  is  zero.  We  shall  then  want  to  find  the  sources  which  put  energy  into 

a  prescribed  radiation  pattern. 

Let  us  now  assume  that  there  are  no  sources  for  t  <  t-.  At  t  "  t^ 

the  sources  are  switched  on,  permitted  to  vary  in  time  in  any  desired  fashion, 

and  then  switched  off  again  at  t  >t-,.  From  (U.l?)  it  is  clear  that  for  t  <  t^ 

we  have  a  solution  of  Maxwell's  equations  without  sources  which  is  a  purely 

transverse  field  which  we  may  write  as 

r         -iep(t-t  ) 
(5.1)  T(x,-t)  »  Y^(xjt)  -  T"   nL(x|p,e)e       ^  Y(p,e,-t^)dp 

e=*l  J 


J 


G(xjt|xStQ)  T'^(x«;t3)dx'  ,     (*  <  t^) 


In  the  time  interval  ^q  <  *  <  t,,  the  transverse  field  is  given  by 
(U.17)  or  (U.19)  . 

The  longitudinal  field  is  given  by  (U.12). 

For  times  t  >  t^  corresponding  to  the  switching  off  again  of  the  sources 
we  have  another  solution  of  Maxwell's  equations  without  sources  which  is  a 
pure  transverse  field; 

(5.2)  I(x;t)  2  lT(x;t)  -  ^I   f  t(x|p,e)e-^'P^*-H^  T(p,e;t  ) 


J 


G(x;t|x'it^)  T^(x';t^)dx'  ,  t  >  t^ 


where 


-iep(t--t«)  -iept,   f  1  ■!  rnf  I 

(5.3)  I(p,e;t^)  -  e     ^  ^  ^Pje^t^)  +  Uni  e     M   J(p,  e;t  •)e"P''  df,   e-*l  . 

^0 
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It  will  be  useful  to  introduce  the  notation 


rh 


(5.U)  F(?,e)  '   J- 


/^  1  r 


t 


Equation  (^.3)  can  be  written 
-iepCt^-t^) 


(6-±l)  . 

-iept^ 


(5.5)  ?(p,e;t, )  -  e     ^-  ^  I(p,6jt^)  ♦  2(2nr/''i  e     ^  F(p,^J  ,   (e-tl)  . 

In  terms  of  the  p-representation,  the  initial  electromagnetic  field  is 

^»* 

given  by  T(p,  e;t^),  while  the  final  one  is  givyn  by  ?(p,  e;t-).  If  we  are 
given  the  sources  J(x;t),  we  can  easily  calculate  the  i  inal  field  fro.ti  the 
initial  field  and  sources  by  using  (5«5)  and  (5»U). 

Let  us  now  consider  the  inverse  problem.  In  this  case  we  are  given 
the  initial  and  final  fields,  Y(p,e;t-.)  and  TCp.fijt^)  respectively,  and  we  are 
required  to  find  the  sources  J(xjt).  We  shall  show  that  the  solution  is  not 
unique,  but  that  we  can  obtain  essentially  unique  results  by  imposing  additional 
conditions. 

From  (5»5)  we  can  find  F(p, e)  from  the  initial  and  final  fields: 
iept,  lept. 


(5.6)  F(p,e) 


2(2n) 


77? 


e     f(p,e;t^)  -e     Y(p,ejtQ) 


where  6*+l  • 

If  we  can  find  ^(p, e^t)  and  hence  J(xjt)  from  F(p, e),  our  problem  is 
solved. 

Let  us  define  the  function  F(p, sjk)  as  being  the  Fourier  transform 
of  J(p, e;t)  with  respect  to  time,  for  e=il. 

That  is 
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(^.7) 


F(p,  e.;k)  >= 


t;e   dt 


< 


0 

,  +00 


J(p,e;t)  -  -^ 


F(p,e;k)e  "^  dk 


-00 


The  condition  that  ^(x;t)  be  real  leads  to  equations  (2,19)  and  hence  also 

(5.7a)  F*(p,6;k)  -  F(  -  p,e;  -  k)  . 

It  is  clear  that  if  we  are  given  the  initial  and  final  fields  and 
hence  F(p, e),  ws  can  take  any  function  F(p, e;k)  which  satisfies  (5«7a)  such  that 

($.8)  F(p,e;ep)  «  F(p,e)  , 

and  obtain  J(p, s;t)  for  e-il  using  the  second  of  equations   (5»7). 
We  can  fiad  a  suitable  real  source  ^(xjt)  frcrni 

(5.9)   J(xjt)  ^Y2.\  I(^!p'^)  I(E.e;^^)dp  » 

irtiere  J(p,+ljt)  is  given  by  (5.7),  f(p»^;t)  is  any  arbitrary  function  which  satisfies 
(2.19)   and  vriiich  vanishes  identically  when  t  <  t-  or  t  >  t,,   and  J(p,0;t)   » 

i  dl  (£,  t-;t) 

p     dt      • 

We  note  that  there  are  two  types  of  lack  of  uniqueness  in  J(x;t)  • 
One  type  consists  in  the  arbitrary  choice  of  function  J(p,  ^'jt)  which 
satisfies  (2.19),  namely  ^(p,  Z^;t)  »  J(-  p,  ^;t),  but  vanishes  identically 
for  times  t  <  t^  and  for  times  t  >  t- .  This  arbitrariness  corresponds  to  the 
possibility  of  having  any  arbitrary  charge  distribution  p(xjt)  for  t-  <  t  <  t^ 
which  would  give  rise  to  the  longitudinal  field  given  by  (1|#12) 
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3x' 


p(x'jt) 

dx« 


x-x' 


When  this  charge  distribution  is  switched  off,  it  will  not  affect  the 
final  transverse  field  in  any  way  whatever.  Thus  the  final  radiation  field 
is  given  only  by  the  sources 


|(x;t)  -21 
c«±l 


'X(x|p,e)  |(£,e;t)dig 

only,  '^o  this  source  we  raay  add  J    where 

|^^(x;t)  =  |l(x!p;t)  |(p,T-,t)dp  +  f%(x|p,0)  $(g,0;t)dg 
which,  when  expressed  in  terms  of  the  arbitrary  charge  density  p  is 
|^^^(x;t)  -  p(x,-t) 


(5.10) 


r^^*) 


8   1 

TT  in? 

3x 


1     9p(x';t) 

J  ~T-S^ ^'  '   (i-1,2,3)  . 


3  af  ^^ 


It  is  clear  that  J^  satisfies  the  equation  of  continuity  £"  r—  "   0, 

2  ,   n  t^o  ax^ 

since  V  f-  —  "  -  6(x). 
^   |x| 

The  choice  of  p(xjt)  and  hence  also  ^   is  not  an  essential  lack 
of  uniqueness  in  the  inverse  problem  because  this  source  does  not  give  rise  to 
transverse  fields  at  any  time:  these  sources  give  rise  only  to  the  longitudinal 
fields  which  have  no  propagation  characteristics. 

The  second  lack  of  uniqueness  in  the  inverse  problem  is  far  more 
important.  It  concerns  the  choice  of  function  F(p, 6;k)  for  e=±l.  Except  for 
the  requirement  that  this  function  satisfy  (5.7a)  and  (5*3),  it  is  arbitrary. 
We  shall  therefore  impose  additional  conditions  which  will  make  the  problem 
unique. 
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6.   The  statement  of  the  inverse  problem  which  leads  to  unique  solutions*  Examples* 
It  will  be  convenient  to  choose  the  origin  of  time  so  that  t^  »  -  T 


arri  t,  «  T  where  T  >  0.  This  can  always  be  done  by  taking  T  -  y(t^-t_)  and 
introducing  a  new  time  coordinate  t'  *  t  -  — = —  .  We  shall  assume  the  new  time 
coordinate  is  always  used  in  what  follows  and  drop  the  prime* 

We  shall  prescribe  the  initial  and  final  fields  in  the  p-representation 
which  are  now  T(p, e;  -T)  and  f(p,e;T), respectively.  Henceforth,  whenever  e 
appears  it  will  be  restricted  to  the  values  6=11. 

The  Ansatz  which  will  lead  to  unique  sources  for  any  choice  of  initial 
and  final  fields  is  the  requirement  that  the  source  vector  J(x;t)  be  represented 
by 

(6.1)   J(xjt)  -  fix)   h^(t)  +  J(x)  h^(t)  , 

where  h  (t)  is  a  prescribed  real  even  function  of  t  and  h  (t)  is  a  prescribed 

real  odd  function  of  t  (the  superscript  U  stands  for  'uneven'  and  is  used 

E        U 
instead  of  0  to  prevent  confusion)*  Neither  h  (t)  nor  h  (t)  is  allowed  to  be 

identically  Jsero.  Our  statement  of  the  inverse  problem  is  that  we  shall  be 

able  to  find  unique  real  U-component  column  vectors  %   (x)  and  ^  (x)  which  are 

functions  of  x  only  for  any  initial  and  final  fields. 

It  will  also  be  shown  that  the  sources  J(3y't)  so  obtained  will  be 

E        U 
independent  of  the  normalization  of  h  (t)  and  h  (t)  in  the  sense  that  if  these 

E         U 
functions  replaced  by  Ah  (t)  and  Bh  (t)  respectively,  where  A  and  B  are  any  real, 

non-zero  constants,  the  same  sources  $(xjt)  will  result. 

£        U 
Let  us  first  introduce  the  Fourier  transforms  of  h  (t)  and  h  (t): 

+T 


(6.2) 


gE(k)  -  -i.     h^(t)e^'^dt  , 


J^k)  .  -L  [   h^t)ei'^dt  . 


-T 
,+T 
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It  is  easily  seen  that  g  (k)  is  a  real  even  function  of  k  and  that 


U 


g  (k)  is  a  purely  Imaginary  odd  function  of  k. 


(6.3) 


g^(k)*  -  g^k) 

g^(-k)  -  g^(k) 

g  (k)  =  -g  (k) 

g"(-k)  -  -g^(k)  . 


From  (?.?)  and  the  relation 


U/        N^U, 


(6.U)   J(g;s;t)  -  |^(p,e)h^(t)  ♦  5^(p,c)h^(t) 


we  see  that 

(6.5)   F(p,6,-k)  -  |%,e)g^(k)  +  J"(p,e)g"(k)  , 

where 


rtl/   s  U/ 


(6.6) 


< 


$^'"(x)  -  r 


X(x|p,6)  f  *  (p,e)dp  . 

r^    »J  ^        **»        ^ 

■CI  TT 

The  reality  conditions  on  J  (x)  and  J  (x)  lead  to  the  requirement  that 

(6.7)  $^"(-P,e)  -|^^*(p,e)  . 

Equation  (5.B),  together  with  (6.3)  yields,  finally, 

(6.8)  F(p,e)  -  J^(p,e)  g^(p)  +  4^(p,6)  g^(p)  . 

Now  it  will  be  convenient  to  introduce  F^(p, e)  and  F  (p, e),  defined  by 

F^(p,6)  -i[F(p,e)  ♦F(-p,6)*] 

(6.9)  {  ^  ^ 
|F"(^t)  -|[f(p,6)  -  F(-p,e)*J  . 
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It  is  clear  that 

F(p,e)  -  F^(p,e)  +  f"(p,6) 


(5.10) 


^k!/' 


<F^(-^,e)  -F^(p,6)* 
F"(-p,e)  -  -F^(p^e)*  . 
Therefore,  using  (6.3)  and  (6.7),  we  have 

F^(p,e)  -  jVe)  g^(p) 

■ 

f"(p,6)  -  6|^(p,e)  g^(p)  . 

P  TT 

Furthermore,  on  using  (5*6),  we  find  the  following  solutions  for  J^  (p»e)  in 
terms  of  the  initial  and  final  fields: 

(6.12)     |Ve) 


(6.11) 


U(2n) 


\j^^     re^^P^Y(p,  e,T)-e-i^PTY(-p,  e,T)*-e-i^P^T(p,  e,-T) 

*e^^PTY(-p,e;-T)j    . 


(6.12a)     f(p,e)  -  Z)l  ^         rei^P'^I(p,e;T)+e-^^P'^T(p,e,-T)*-e-i'^PTT(^,.,-T) 

U(2n)^/V(p)   L  -- 

-e^^P^(-p,e,Ti*]    . 


Finally,   we  have  the  sources 


(6.13)     f(x;t)  -  r       [    X(x|p,e)  fV^Mp  h^(t)  ♦  ^         X(xlp,e)  fVe)dp  h"(t)  . 


We  see  that  the  presence  of  the  ratios  ■  w^   ,  u.   '  in  (6.13)  when 


E,U, 


g  (p)   g  (p) 


(6.12)  and  (6.12a)  are  substituted  for  J^  (p»e)  leads  us  to  our  conclusion 

E         tJ 
that  J{xjt)  is  independent  of  the  normalizations  of  h  (t)  and  h  (t)  in  the 

sense  mentioned  before. 

To  siimmarize:  we  prescribe  the  initial  field  T(p,  e;-T)  and  the  final 
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field  Y(p,  ejT).  We  require  Kxjt)  to  have  the  form  (6.1)  where  h  (t)  and 

h  (t)  are  given  real  even  and  odd  functions  respectively  of  t.  We  then  compute 

f^(p,e)  and  f"(p,e)  using  (6.12)  and  (6.12a)  and  finally  |(xjt)  fron  (6.13). 

It  should  be  noted  that  it  is  only  for  exceptional  cases  of  initial 
and  final  fields  that  either  f  (x)  or  f^(x)  is  identically  zero  or  that  \   (x)-f^(x). 
These  situations  correspond  to  the  cases  in  which  the  sources  can  be  written 
J^(x)h^(t),  J^(x)h^(t),  or  |(x)h(t)  where  h(t)  is  an  arbitrary  given  function  of  t. 

We  shall  now  consider  two  examples  of  the  inverse  problem. 

Example  1.  We  shall  take  the  field  before  the  sources  are  switched  on 
to  be  zero.  After  the  sources  are  switched  off  the  field  is  to  consist  of  a 
circularly  polarized  wave  travelling  in  the  positive  x-directton  with  frequency  v. 

Hence 


(6.1h) 


T(g,-e,-T)  :  0 


|Y(p,e;T)  -  K£(p^-  v)  6(Py)  5(Pjj)  6^^^^  ,  v  >  0  . 


(6.15) 


Furthermore,  we  shall  take 


E 


h''(t)  -  A5(t) 
h^(t)  -  B6'(t)  • 

E        U 
All  sufficiently  short  time  h  (t)  and  h  (t)  functions  can  be  approximated 

by  the  functions  given  in  (6.15).  Furthermore,  we  may  take  T  to  be  arbitrarily 

small.  In  fact  we  shall  take  T  to  be  zero  after  the  various  integrations  over 

time  have  been  performed. 

Then,  in  terms  of  the  x-representation,  the  field  is  identically  zero 

for  t  <  0  and  for  t  >  0  is  given  by  (see  3,11) 
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(6,36) 


E  -  0  , 


y   >/?  (2n) 


v^  (2n) 


"x 


sin  v(x-t)  , 


^   ^  ■-  cos  V(x-t)  ,        H^ "^  K/x- 

j-r^  sin  v(x-t)  ,       H^  »  -^^— -^575-  cos  v(x-t)  . 


^  '  y?  (2n)^/2 


The  calculations  for  the  sources  are  quite  straightforward  and  one 


obtains 


0 
0 


(6.17)     f  (x)  -  -r-TT^       -  cos  vx   ,     f  (x) 


K 


A  16  n 


sin  vx 


B  l6n 


^v 


0 

0 

sin  vx 

cos  vx. 


and 


(6.38) 


¥V^^  '  ^^ 


-  cosvx  5(tj+  SiiLH  6.(t) 
sinvx  6(t)+  £22_H  6t(t) 


One  can  add  to  this  source  the  arbitrary  charge  distribution  p(xjt) 
which  we  may  choose,  if  we  wish,  to  take  the  form 

p(x;t)  -  p(x)6(t) 
where  p(x)  is  arbitrary.  This  given  rise  to  the  arbitrary  additional  sources 

'p(x)6(t) 


1  d 
Hn  5x 


p(x') 


(6.19) 


*Arb,  .X  1   1  a  1 


1^^'  I 

|x-xM 
P(x') 


dx'6'(t) 


dx'6«(t) 


v^^)t£t"-^''*^ 


Having  obtained  the  sources,  we  can  now  obtain  the  fields  for  all  time 
using  techniques  given  in  Section  U.  The  transverse  field  is  given  by 
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(6.20)   f^(x;t)  -  ^-.5-75. 

y?  (2n)^/2 


^J(t) 


iv(x-t) 


where  7]  (t)  is  the  Heaviside  function 

(6.20a)    y  (t)  -  0   t  <  0 
-  1   t  >  0  . 

The  longitudinal  field  is 
'0 


(6.21)   Y^(3yt) 


-  i 


dx»5(t) 


Example  2 . 

In  the  previous  example  the  final  field  filled  all  of  space  shortly  after 
the  sources  were  switched  off.  As  a  consequence  it  is  seen  that  the  sources 
also  fill  all  of  space  during  the  time  that  they  are  on. 

In  the  present  example  we  shall  study  the  case  in  which  the  time 
dependence  of  the  sources  is  the  same  as  before  and  where  we  again  require 
the  initial  field  to  be  identically  zero.  In  contrast  to  the  previous  example, 
however,  we  shall  require  the  final  transverse  field  to  be  highly  concentrated  in 
space  immediately  after  the  sources  are  switched  off.  We  shall  then  see  that  as 
long  as  the  sources  are  on,  they  will  also  be  highly  concentrated  in  space. 

It  would  be  nice  to  consider  the  field  immediately  after  the  switch-off 

to  be 

(6.22)   T^(x)  =■  R6(x) 

where  R  is  the  f our-ccwiponent  vector 

'0 


(6.22a) 
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vhere  R.  is  any  complex  number  R.  ■  L.+  iM..  However,  I  ,  as  it  stands,  is  not 
suitable  for  a  final  field  because  it  is  not  purely  transverse.  Hence  we  shall 
subtract  the  part  whose  divergence  is  not  zero. 

In  accordance  with  the  previous  sections  we  may  expand  i  (x)  as 


JX(x|p,i 


(6,23)  I^(x)  -r  )M^\Pft)   T*(p,e)dp  +   X(x|p,0)  T^(p,0)dp  . 
Now,  the  part  for  which  the  divergence  does  not  vanish  is 


(6.2U)  T^(x)  =   X(x|p,0)  T^(p,0)dp  . 


Moreover,  since 

(6.25)  T^(P,0)  -  r  J  Xi(xlp,0)  Tj(x)d^ 

we  have,  on  using  the  explicit  form  of^  .(x|p,0), 

1  ~  /v 

(6.26)  Y^^(x)  =  r  [  dx'T^(x')  t.(x|p,0)X*(x'|p,0)dp 

^nce,  we  shall  take  as  the  field  immediately  after  the  sources  are 
switched  off 

(6.27)  I,  -  R,6(x)  ♦  1.  ^  Rj  -ji-  ^  ,   (i.1,2.3)  . 

which  is  still  highly  localized  near  the  origin  and  is  now  a  purely  transverse  wave, 
The  sources  which  give  rise  to  (6.2?)  are  easily  calculated: 


(6.28) 


llM   -  0 


3       2 

f^M   .  4v  M.6(x)  +  4-  r  M,  -^ — r  -rr 


(i-1,2,3) 
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(6.28a) 


r 


(Un)'B  L 


T   9    I   a  "1  1 


^2  n 


■"l^lTt 


w 


■3  -"    (Un)^B  L^  9x" 
In  addition  to  these  sources,  one  can  add  the  arbitrary  sources  given 
by  (6.19). 

For  t  >  0,  it  is  easy  to  show  that  the  field  is  given  by 

(6.29)   Y. (x;t)  -  F  G^j(x;t|0;O)Rj 

where  G  is  the  Green's  function  given  by  (3.17)  or  (3.18). 
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